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This  paper  illustrates  the  value  of  SeaWiFS  ocean  color  imagery  in  assessing 
the  ability  of  three  data-assimilative  ocean  models  (configured  in  five  predic¬ 
tion  systems)  to  map  mesoscale  variability  in  the  Gulf  of  Mexico  (i.e.,  the  Loop 
Current  and  associated  warm  and  cold  eddies)  and  in  helping  to  diagnose  specific 
strengths  and  weaknesses  of  the  systems.  In  addition,  the  study  clearly  illustrates 
that  biological  responses  of  the  surface  waters  are  strongly  linked  to  the  physical 
events  and  processes. 


1.  INTRODUCTION 

An  accurate  depiction  of  currents  and  associated  features 
such  as  mesoscale  eddies,  fronts,  and  jets,  is  essential  to  any 
ocean  forecasting  system.  The  ability  to  use  high-resolu¬ 
tion  ocean  models  for  data  assimilation  and  forecasting  is 
critical  to  the  development  of  these  nowcast/forecast  systems. 
Applications  include  assimilation  and  synthesis  of  global 
satellite  surface  data,  ocean  prediction,  optimum  track  ship 
routing,  search  and  rescue,  antisubmarine  warfare  and  surveil¬ 
lance,  tactical  planning,  high  resolution  boundary  conditions 
that  are  essential  for  even  higher  resolution  coastal  models, 
sea  surface  temperature  for  long  range  weather  prediction, 
inputs  to  shipboard  environmental  products,  environmental 
simulation  and  synthetic  environments,  observing  system 
simulations,  ocean  research,  pollution  and  tracer  tracking, 
inputs  to  water  quality  assessment,  and  so  forth. 

Since  most  high  horizontal  resolution  observational  data 
(i.e.,  altimetry  and  sea  surface  temperature)  are  assimilated 
in  these  systems,  there  are  very  few  independent  observa- 
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tions  with  sufficient  coverage  to  be  useful  in  assessing  the 
ability  of  ocean  prediction  models  to  map  mesoscale  ocean 
features.  In  this  paper,  we  assess  the  performance  of  several 
data  assimilative  ocean  model  systems  against  indepen¬ 
dent  ocean  color  measurements  in  the  Gulf  of  Mexico.  The 
near-surface  biological  and  optical  properties,  derived  from 
ocean  color  imagery,  provide  an  excellent  evaluation  method 
for  assessing  the  capability  of  data-assimilative  numerical 
models  to  define  the  spatial  and  temporal  mesoscale  vari¬ 
ability.  Although  these  properties  are  non-conservative 
and  are  influenced  by  biological  processes,  these  optical 
gradients  are  an  independent  data  set  with  the  ability  to 
track  and  monitor  mesoscale  features  such  as  warm  and  cold 
core  rings,  coastal  filaments  and  the  Loop  Current  position. 
Feature  tracking  using  a  time  series  of  these  independent 
observations  in  the  Gulf  of  Mexico  is  a  very  effective  way 
to  compare  different  ocean  prediction  systems. 

Because  of  the  logistics  associated  with  routine  analysis 
and  prediction,  most  large-scale  ocean  prediction  systems 
are  located  within  government  centers.  Recently,  the  Global 
Ocean  Data  Assimilation  Experiment  (GODAE)  has  pro¬ 
vided  an  international  framework  for  such  systems  to  be 
evaluated  and  inter-compared.  The  rationale  for  restricting 
ourselves  to  the  systems  compared  in  this  paper  (see  section 
3  for  details)  is  that  they  all  use  the  same  atmospheric  forc¬ 
ing  and  assimilate  the  same  observational  data.  This  allows 
us  to  evaluate  the  performance  of  the  systems  independently 
of  the  choice  made  for  external  inputs  (atmospheric  forcing 


20060619008 


88  ASSESSMENT  OF  DATA  ASSIMILATIVE  OCEAN  MODELS  IN  THE  GULF  OF  MEXICO  USING  OCEAN  COLOR 


or  ocean  data  input).  Another  example  of  a  Gulf  of  Mexico 
data  assimilative  system  is  discussed  by  Kantha  et  al.  in 
this  volume. 

2.  OCEAN  COLOR 

The  data  set  used  to  evaluate  the  forecasting  systems  con¬ 
sists  of  satellite  measurements  of  ocean  color  data  obtained 
by  the  Sea-viewing  Wide  Field-of-view  Sensor  (SeaWiFS)  in 
orbit  on  the  OrbView-2  (formerly  SeaStar)  platform.  Ocean 
color  satellites  such  as  SeaWiFS  and  MODIS  provide  new 
capability  to  quantitatively  monitor  biological  and  optical 
properties  of  the  near  surface  ocean  [Arnone  and  Parsons. 
2004].  Accurate  sensor  calibration  and  atmospheric  cor¬ 
rection  are  required  in  order  to  detect  the  subtle  changes  in 
ocean  color  which  are  used  to  uncouple  the  in-water  com¬ 
ponents. 

The  apparent  color  of  the  ocean,  which  can  range  from 
deep  blue  to  varying  shades  of  green  and  ruddy  brown,  is 
primarily  determined  by  the  concentration  of  substances 
and  particles  in  the  euphotic  (lighted)  zone  of  the  upper 
ocean.  Living  phytoplankton  (which  contain  chlorophyll  and 
associated  photosynthetic  pigments),  inorganic  sediments, 
detritus  (particulate  organic  matter),  and  dissolved  organic 
matter  all  contribute  to  the  color  of  the  ocean.  When  visible 
light  from  the  sun  illuminates  the  ocean  surface,  it  is  subject 
to  several  optical  effects.  Foremost  among  these  effects  are 
light  absorption  and  scattering.  Spectral  absorption  selec¬ 
tively  removes  some  wavelengths  of  light  while  allowing 
the  transmission  of  other  wavelengths.  In  open  ocean  waters 
absorption  is  primarily  due  to  the  photosynthetic  pigments 
(chlorophyll)  present  in  phytoplankton,  whereas  in  coastal 
waters  additional  absorption  processes  occur  from  the 
absorption  of  detritus  and  colored  dissolved  organic  matter 
(CDOM).  Following  spectral  absorption,  backscattercd  light 
(resulting  from  organic  and/or  inorganic  particulate  matter 
suspended  in  the  water)  produces  a  modified  light  radiating 
from  the  ocean  surface,  the  “water-leaving  radiance."  i.e..  the 
ocean  color.  Spacebornc  radiometers  such  as  SeaWiFS  mea¬ 
sure  the  spectral  radiance  intensity  at  specific  wavelengths. 
These  measurements  are  used  to  determine  quantitative 
constituents  in  the  water  column  that  interact  with  visible 
light,  such  as  chlorophyll.  Additional  components  of  the 
absorption  and  backscattering  light  field  can  be  uncoupled 
from  the  ocean  color  signatures,  but  are  beyond  the  scope  of 
this  paper.  We  used  the  chlorophyll  component  in  this  study 
because  it  is  the  dominant  contributor  to  the  ocean  color  in 
open  oceamprocesses. 

The  chlorophyll  products  from  SeaWiFS  represent  the  inte¬ 
grated  near-surface  concentrations  within  the  first  attenu¬ 
ation  length  of  visible  light  (490  nm).  For  Loop  Current 


waters,  this  is  approximately  20  meters,  whereas  in  turbid 
coastal  areas,  this  is  reduced  to  as  little  as  two  to  three 
meters.  This  product  is  different  from  a  satellite  SST  product 
that  only  senses  the  skin  temperature  and  responds  to  diurnal 
heating.  Ocean  color  products  can  be  used  to  detect  ocean 
features  throughout  the  year,  which  is  especially  useful  dur¬ 
ing  the  summer  months  when  the  SST  in  the  Gulf  of  Mexico 
is  nearly  isothermal  and  cannot  be  used  for  feature  detection. 
However,  chlorophyll  concentrations  can  change  rapidly 
through  phytoplankton  growth  and  decay  and  through  verti¬ 
cal  movements  of  phytoplankton  into  or  out  of  view  of  the 
satellite. 

The  SeaWiFS  mission  is  a  public-private  partnership 
between  NASA  and  Orbital  Sciences  Corporation.  SeaWiFS 
was  launched  in  August  1997  and  acquires  approximately 
15  pole-to-pole  orbital  swaths  (~  2800  km  wide)  of  data 
per  day,  and  approximately  90%  of  the  ocean  surface  is 
scanned  every  two  days.  The  SeaWiFS  ocean  data  set  used 
in  this  effort  was  collected  and  processed  daily  for  the  Gulf 
of  Mexico  at  the  Naval  Research  Laboratory  ( NRL)  receiv¬ 
ing  site  using  the  Automatic  Processing  System  developed 
for  the  Naval  Oceanographic  Office  to  prov  ide  operational 
products.  Atmospheric  correction  and  in-water  chlorophyll 
algorithms  (OC4)  were  applied  [Arnone  el  al..  1998;  O'Reilly 
et  al..  1998]  as  extensions  of  the  NASA  algorithms  to  address 
the  coastal  waters.  Note  that  the  OC4  algorithm  uses  a  simple 
ratio  of  spectral  channels  at  443.  488.  520  and  550  nm  of 
the  water-leaving  radiance  to  determine  chlorophyll.  This 
simple  algorithm  is  closely  associated  with  total  absorption 
coefficient,  which  in  the  open  ocean  is  closely  associated 
with  chlorophyll  concentrations  [Lever  at..  1998],  However, 
for  other  locations,  such  as  in  coastal  areas,  the  algorithm 
significantly  overestimates  chlorophyll,  and  the  values  are 
closely  associated  with  the  total  absorption  coefficient  rather 
than  chlorophyll. 

Because  the  majority  of  visible  imagery  is  limited  by  cloud 
cover,  the  latest  pixel  composite  is  created  daily  to  represent 
the  optimum  coverage  of  ocean  features  within  the  last  ~ 
days.  The  compositing  technique  provides  the  most  recent 
image  of  ocean  conditions  each  day.  A  seven-day  period 
effectively  removes  cloud  cover  and  still  retains  consistency 
in  feature  tracking  to  monitor  the  slower  mesoscale  ocean 
features  for  comparison  with  the  model  results. 

Ocean  color  can  change  rapidly  in  response  to  I)  open 
ocean  biological  processes  occurring  at  daily  time  scales 
and  2)  subduction  of  surface  chlorophyll  below  the  satellite 
detection  depth.  Phytoplankton  concentrations  can  double 
within  a  day  provided  optimum  light  and  nutrient  avail¬ 
ability.  More  traditionally,  chlorophyll  response  in  upvvell- 
ing  regions  occurs  on  time  scales  of  4-7  days  as  the  bloom 
matures,  and  blooms  with  a  continuous  nutrient  supply  can 
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remain  observable  for  months.  A  similar  response  occurs  in 
frontal  boundaries  where  a  steady  increase  in  chlorophyll  is 
observed  with  constant  frontal  upwelling.  Similarly,  blooms 
can  degrade  rapidly  in  4-7  days  as  nutrients  are  rapidly 
depleted.  However,  these  daily  color  changes  mainly  influ¬ 
ence  the  magnitude  of  the  chlorophyll  concentration,  whereas 
the  location  of  the  fronts  and  eddies  is  driven  primarily  by 
the  circulation.  Therefore,  the  use  of  the  seven  day  compos¬ 
ite  provides  a  reasonable  location  of  the  mesoscale  features 
which  are  identified  by  the  models. 

3.  THE  DATA  ASSIMILATIVE  OCEAN  MODELS 

The  three  data  assimilative  ocean  models  (configured  into 
five  prediction  systems)  discussed  in  this  paper  are  designed 
for  nowcasting  and  forecasting  of  the  mesoscale,  large  scale, 
and  upper  ocean.  Two  of  the  prediction  systems  use  the  NRL 
Layered  Ocean  Model  (NLOM)  [Hurlburt  and  Thompson , 
1980;  Wallcraft  et  al.,  2003]  and  two  use  the  Navy  Coastal 
Ocean  Model  (NCOM)  [Barron  et  al.,  2004, 2006,  Kara  et  al, 
2005],  all  four  developed  at  the  Naval  Research  Laboratory, 
Stennis  Space  Center.  The  final  system  is  based  on  HYCOM 
(H  Ybrid  Coordinate  Ocean  Model)  [Bleck,  2002;  Chassignet 
et  al.,  2003;  Halliwell,  2004]  which  has  been  developed  by 
a  community  effort  (http://hycom.rsmas.miami.edu)  sup¬ 
ported  by  the  National  Ocean  Partnership  Program  (NOPP). 
The  rationale  for  restricting  ourselves  to  these  systems  is 
that  they  all  use  the  same  atmospheric  forcing  from  the 
Fleet  Numerical  Meteorology  and  Oceanography  Center 
(FNMOC)  Navy  Operational  Global  Atmospheric  Prediction 
System  (NOGAPS)  [Hogan  and  Rosmond.  1991;  Rosmond 
et  al.,  2002],  and  assimilate  the  same  observational  data,  the 
altimeter  sea  surface  .height  (SSH)  as  well  as  the  operational 
1/8°  MODAS  [Modular  Ocean  Data  Assimilation  System; 
Fox  et  al.,  2002}  sea  surface  temperature  (SST)  analyses  of 
AVHRR  (Advanced  Very  High  Resolution  Radiometer)  data. 
This  allows  us  to  evaluate  the  performance  of  the  systems 
independently  of  the  choice  made  for  external  inputs  (atmo¬ 
spheric  forcing  or  ocean  data  input). 

3. 1  The  Numerical  Models 

NLOM  was  designed  for  use  in  a  first  generation  global 
ocean  prediction  system  by  allowing  high  horizontal 
resolution  at  much  lower  computational  cost  than  other 
models  [Hurlburt,  1984],  NLOM  has  been  configured 
globally  from  72°S  to  65°N  at  1/16°  and  1/32°  (~7  and  3.5  km 
mid-latitude,  respectively)  resolution.  However,  NLOM  has 
only  7  layers  in  the  vertical,  including  the  mixed  layer,  and 
excludes  the  Arctic  and  most  shallow  water  (i.e.,  the  lateral 
boundaries  usually  follow  the  200  m  isobath.). 


Global  NCOM  was  designed  to  complement  NLOM  by 
covering  these  regions  and  using  much  higher  vertical 
resolution,  but  with  coarser  1/8°  (15-16  km  mid-latitude) 
horizontal  resolution.  NCOM  has  a  mixed  a- z  vertical 
coordinate  [Barron  et  al.,  2005]  where  terrain-following 
a  coordinates  are  used  when  the  bottom  depth  is  <137  m 
in  the  current  global  configuration.  This  configuration 
has  40  levels  in  the  vertical  with  vertical  resolution  ~1  m 
very  near  the  surface  for  fine  resolution  of  the  mixed  layer. 
Global  NCOM  encompasses  the  open  ocean  to  5  m  depth 
in  a  curvilinear  global  model  grid  with  1/8°  latitudinal  grid 
spacing  at  45°N.  The  grid  extends  from  80°S  to  a  complete 
Arctic  cap  with  grid  singularities  mapped  into  Canada  and 
Russia.  A  regional  configuration  of  NCOM  is  discussed  in 
section  5.  . 

Recent  model  comparison  exercises  performed  in  Europe 
(DYnamics  of  North  Atlantic  MOdels — DYNAMO) 
[Willebrand  et  ai,  2001]  and  in  the  U.S.  (Data  Assimilation 
and  Model  Evaluation  Experiment — DAMEE)  [Chassignet 
et  al.,  2000]  have  shown  that  no  single  vertical  coordinate 
(depth,  density,  or  terrain-following  sigma)  can  by  itself 
be  optimal  everywhere  in  the  ocean.  Isopycnal  (density 
tracking)  layers  are  best  in  the  deep  stratified  ocean,  z-Ievels 
(constant  fixed  depths)  are  best  used  to  provide  high  vertical 
resolution  near  the  surface  within  the  mixed  layer,  and 
sigma-levels  (terrain-following)  are  often  the  best  choice 
in  shallow  coastal  regions.  HYCOM  combines  all  three 
approaches  by  choosing  the  optimal  distribution  at  every 
time  step.  The  hybrid  coordinate  is  one  that  is  isopycnal 
in  the  open,  stratified  ocean,  but  that  makes  a  dynamically 
smooth  transition  (via  the  layered  continuity  equation)  to 
terrain-following  coordinates  in  shallow  coastal  regions, 
and  to  pressure  coordinates  in  the  mixed  layer  and/or 
unstratified  seas.  The  capability  of  assigning  additional 
coordinate  surfaces  to  the  oceanic  mixed  layer  in  HYCOM 
provides  the  option  of  implementing  sophisticated  vertical 
mixing  turbulence  closure  schemes  (see  Halliwell  [2004]  for 
a  review).  HYCOM  has  been  configured  globally,  on  basin 
scales,  and  regionally. 

The  1/12°  (~7  km  mid-latitude  resolution)  North  Atlantic 
HYCOM  system  is  configured  from  28°S  to  70°N,  including 
the  Mediterranean  Sea  (see  Chassignet  et  al.  [2005]  for 
details).The  vertical  resolution  consists  of  26  hybrid  layers, 
with  the  top  layer  typically  at  its  minimum  thickness  of  3 
m  (i.e.,  in  fixed  coordinate  mode  to  provide  near-surface 
values).  In  coastal  waters,  there  are  15  sigma-levels,  and 
the  coastline  is  at  the  10  m  isobath.  The  northern  and 
southern  boundaries  are  treated  as  closed,  but  are  outfitted 
with  3°  buffer  zones  in  which  temperature,  salinity  and 
pressure  are  linearly  relaxed  toward  their  seasonally  varying 
climatological  values. 
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3.2  Data  Assimilation 

The  1/16°  and  1/32°  global  NLOM-based  systems  [Smedstad 
et  al.,  2003;  Shriver  et  al,  2005]  assimilated  track  data  from 
available  real-time  satellite  altimeter  observations  (GFO 
and  JASON)  using  an  optimal  interpolation  (OI)  deviation 
SSH  analysis  with  the  model  field  as  the  first  guess  and 
mesoscale  covariances  calculated  from  Topex-Poseidon 
and  ERS-2  altimeter  data  [Jacobs  et  al.,  2001].  A  statistical 
inference  technique  [Hurlburt  et  al.,  1990]  updates  all  layers 
of  the  model  based  on  the  analysed  SSH  deviations,  including 
geostrophic  updates  of  the  velocity  field  outside  an  equatorial 
band.  The  global  model  is  then  updated  to  produce  a  nowcast 
using  slow  insertion  to  further  reduce  gravity  wave  generation, 
fit  addition,  model  deviations  from  full  field  MODAS  SST 
analyses  are  assimilated  in  the  form  of  heat  fluxes. 

The  1/8°  global  NCOM  system  [Barron  et  al.,  2004.  2005] 
uses  relaxation  (i.e.,  nudging)  to  assimilate  3-D  synthetic 
temperature  and  salinity  fields  generated  from  MODAS  SST 
analyses  and  the  1/16°  NLOM  SSH  fields  using  the  MODAS 
system  algorithms  and  statistics  [Fox  et  al.,  2002],  The 
nudging  in  NCOM  is  weak  from  below  the  surface  to  the 
base  of  the  mixed  layer  to  allow  greater  model  impact  on  the 
mixed  layer. 

The  data  assimilation  in  the  current  1/12°  North  Atlantic 
HYCOM  system  consists  of  assimilating  daily  operational 
MODAS  SSH  analyses  of  available  real-time  satellite  altim¬ 
eter  observations  (GFO  and  JASON)  [Chassignet  et  al.. 
2005].  The  Cooper  and  Haines  [1996]  technique  is  used  to 
project  the  surface  information  to  the  interior  of  the  ocean.  A 
relaxation  to  the  MODAS  SST  analysis  is  also  included. 

All  systems  provide  nowcasts  and  forecasts  at  least  once 
,  a  week,  with  forecasts  for  30  days  from  the  global  NLOM 
systems,  5  days  from  the  global  NCOM  system,  and  14  days 
from  the  1/12°  Atlantic  HYCOM  system. 

3.3  The  Mean  Sea  Surface  Height 

In  order  to  assimilate  the  anomalies  determined  from 
the  satellite  altimeter  data  into  the  numerical  model,  it  is 
necessary  to  know  the  oceanic  mean  SSH  over  the  time 
period  of  the  altimeter  observations.  Unfortunately,  the  geoid 
is  not  known  accurately  on  the  mesoscale.  Several  satellite 
missions  are  either  underway  or  planned  to  try  to  determine 
a  more  accurate  geoid,  but  until  the  measurements  become 
accurate  to  within  a  few  centimeters  on  scales  down  to 
approximately  30  km.  one  has  to  define  a  mean  oceanic  SSH. 
At  the  scales  of  interest  (tens  of  kilometres),  it  is  necessary 
to  have  the  mean  of  major  ocean  currents  and  associated 
SSH  fronts  sharply  defined.  This  is  not  feasible  from  coarse 
hydrographic  climatologies  (~l“  horizontal  resolution)  and 


the  approach  taken  by  the  above  systems  has  been  to  use  a 
model  mean  SSH.  This  requires  a  fully  eddy-resolving  ocean 
model  which  is  consistent  with  hydrographic  climatologies 
and  with  fronts  in  the  correct  position. 

For  the  NLOM-based  system,  the  mean  SSH  from 
a  global  NLOM  1/16°  run  without  data  assimilation  was 
compared  to  the  mean  dynamic  height  calculated  from 
available  temperature  and  salinity  measurements  and  then 
modified  with  data  over  the  period  of  the  satellites.  To  help 
determine  the  position  of  the  fronts,  mean  frontal  positions 
from  satellite  IR  plus  the  variability  from  satellite  altimeter 
data  were  also  used.  Following  this  effort  to  accurately 
determine  the  mean  position  of  the  front,  the  model  mean 
was  modified  by  moving  the  SSH  features  in  an  elastic  way 
(rubber-sheeting)  [Hord,  1982;  Clarke,  1990;  Carnes  at  al., 
1996]  (see  Smedstad  et  al.  [2003]  for  an  illustration).  The 
1/8°  global  NCOM  uses  the  same  mean  SSH  and  assimilates 
the  1/16°  NLOM  SSH  fields.  The  HYCOM-based  forecasting 
system  uses  a  different  model  mean,  one  generated  by  a 
previous  1/12“  North  Atlantic  simulation  performed  with  the 
Miami  Isopycnic  Coordinate  Ocean  Model  (MICOM)  (see 
Chassignet  and  Garraffo  [2001]  for  a  discussion). 

4.  COMPARISON  TO  OCEAN  COLOR 

Plates  1-4  compare  SSH  analyses  from  the  three  ocean 
data  assimilative  systems  which  overlie  ocean  color  imagery 
from  SeaWiFS  in  the  Gulf  of  Mexico.  Each  Plate  is  fora  dif¬ 
ferent  date  (July  9.  July  19.  July  26.  and  August  8.  2003)  and 
contains  results  from  the  I  lh°  global  NLOM  (8  km  resolu¬ 
tion  in  the  Gulf  of  Mexico),  the  I  32°.gIobal  NLOM  (4  km 
resolution),  the  I  8°  global  NCOM  (17  km  resolution),  and 
the  1, 12°  Atlantic  HYCOM  (8  km  resolution).  These  dates 
correspond  to  a  Loop  Current  shedding  event  as  well  as 
further  interactions  between  the  newly  formed  ring  and  the 
Loop  Current.  These  snapshots  clearly  illustrate  how  ocean 
color  has  the  potential  to  characterize  the  performance  of 
each  system,  but  do  not  allow  for  a  quantitative  assessmen: 
of  the  systems.  The  SeaWiFS  imagery  is  a  composite  of 
the  most  recent  cloud-free  pixel  over  the  preceding  6  days 
and  is  colored  prismatically  from  low  (violet)  to  high  (red). 
Both  bright  and  dark  areas  of  ocean  color  provide  useful 
information.  Dark  areas  of  low  chlorophyll  (<  .1  nig/m3) 
usually  correspond  to  Caribbean  water  from  outside  the 
Gulf  of  Mexico  and  can  be  seen  inside  the  Loop  Current 
and  shed  eddies.  When  the  model  is  in  good  agreement  with 
SeaWiFS  imagery,  the  SSH  contours  of  the  Loop  Current 
and  shed  eddies  capture  the  dark  areas  quite  well.  The  mis¬ 
matches  show  up  clearly  and  the  four-way  comparison  is  very 
effective  in  bringing  them  out.  One  has  to  be  careful  when 
performing  the  evaluation,  since  the  overlain  SSH  contours 
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Plate  I.  SSH  from  four  ocean  data  assimilative  systems  on  09  July  2003  (contour  interval  of  10  cm  and  normalized 
over  the  displayed  area),  overlying  ocean  color  imagery  from  SeaWiFS  (natural  log  of  the  concentration  in  mg  mJ). 

The  model  land-sea  boundary  is  defined  by  the  black  contour  line. 
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Ocean  Color  vs.  1/16°  NLOM  SSH  - 19  July  2003 
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Plate  2.  SSH  from  four  ocean  data  assimilative  systems  on  19  July  2003  (contour  interval  of  10  cm  and  normalized 
over  the  displayed  area),  overlying  ocean  color  imagery  from  SeaVViFS  (natural  log  of  the  concentration  in  mg  mJ). 
The  model  land-sea  boundary  is  defined  by  the  black  contour  line. 
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can  occasionally  obscure  some  of  the  features  seen  in  the 
imagery.  The  bright  areas  of  high  chlorophyll  (>.7mg/m3), 
usually  generated  along  the  coast  in  the  Gulf  of  Mexico,  tend 
to  be  advected  into  plumes  by  strong  currents  but  sometimes 
are  also  present  in  the  center  of  cyclonic  (counterclockwise) 
eddies. 

09  July  2003:  Plate  1  clearly  illustrates  the  final  process 
leading  to  the  separation  of  a  Loop  Current  ring.  Specifi¬ 
cally,  there  is  a  bright  patch  of  high  chlorophyll  (~.5mg/m3) 
that  is  being  advected  across  the  Loop  Current  as  the  eddy 
separates.  The  four  models  differ  from  each  other  signifi¬ 
cantly.  The  NLOM  1/16°  does  a  very  good  job  of  capturing 
the  complex  pattern  of  the  Loop  Current,  especially  south¬ 
east  of  the  bright  chlorophyll  patch.  The  soon-to-be  shed 
eddy,  however,  does  not  extend  as  far  north  as  the  SeaWiFS 
imagery  suggest.  The  NLOM  1/32°  does  a  better  job  in  that 
regard,  but  the  westward  front  of  the  Loop  Current  does  not 
extend  sufficiently  southwest.  The  wedge  represented  by  the 
bright  chlorophyll  patch  is  also  not  well  defined  in  NLOM 
1/32°.  Close  examination  of  NLOM  1/16°  and  NLOM  1/32° 
also  show  that  the  westward  extent  of  the  eastern  front  of 
the  Loop  Current  closest  to  the  Florida  Keys  may  be  slightly 
exaggerated  in  NLOM  1/16°.  Despite  the  fact  that  NCOM 
1/8°  assimilates  the  SSH  fields  of  NLOM  1/16°,  the  two 
nowcasts  differ  significantly.  NCOM  1/8°  clearly  misses  a 
large  portion  of  the  northward  extent  of  the  Loop  Current  and 
shows  a  Loop  Current  ring  already  well  defined.  Part  of  the 
bright  patch  of  chlorophyll  can  actually  be  seen  inside  the 
newly  formed  ring,  a  clear  indication  of  the  misplacement 
of  the  frontal  features.  The  wedge  represented  by  the  bright 
chlorophyll  patch  is  perhaps  best  represented  by  HYCOM 
1/12°  but,  as  in  NLOM  1/32°,  the  western  front  of  the  Loop 
Current  does  not  extend  sufficiently  far  south. 

19  July  2003:  This  snapshot  (Plate  2)  is  especially  illustrative, 
since  the  SeaWiFS  imagery  clearly  shows  a  separated  Loop  Cur¬ 
rent  ring.  Furthermore,  the  irregular  shape  of  the  ring  is  quite  chal¬ 
lenging  for  the  models.  For  that  specific  day,  NLOM  1/32“  agrees 
best  with  the  color  measurements.  Of  particular  note  is  the  accurate 
representation  of  the  shape  and  position  of  the  ring  and  of  the  Loop 
Current  as  depicted  by  the  low  chlorophyll  water  (deep  purple).  The 
HYCOM  1/12°  also  does  quite  well  at  representing  these  features, 
but  the  ring  does  not  extend  far  enough  to  the  southwest  and  the 
Loop  Current  does  not  extend  far  enough  to  the  north.  The  Loop 
Current  is  pretty  well  represented  in  both  NLOM  1/16°  and  NCOM 
1/8°  (except  at  the  neck  in  NLOM  1/16°),  but  both  systems  fail  to 
place  the  newly  formed  ring  properly.  The  location  of  the  circular 
ring  in  NCOM  1/8"  (already  formed  on  July  09, 2003)  shows  a  large 
mismatch  between  the  model  results  and  the  chlorophyll  observa¬ 
tions.  Additionally,  only  NLOM  1/32“  and  HYCOM  1/12"  capture 


the  small  cyclonic  eddy  depicted  by  high  chlorophyll  on  the  south 
side  of  the  Loop  Current  ring.  HYCOM  1/12°  is  the  most  success¬ 
ful  in  depicting  a  large  anticyclonic  eddy  in  the  western  Gulf  with 
NCOM  1/8°  in  second  place.  NLOM  1/32°  is  the  least  successful. 

26  July  2003:  The  Loop  Current  ring  that  showed  as  clearly 
separated  on  July  19,  2003  is  still  separated,  but  showing 
signs  of  being  recaptured  by  the  Loop  Current  (Plate  3). 
NCOM  1/8°  shows  very  little  change  from  09  July  2003  and 
fails  completely  to  move  the  ring  northeastward  (actually 
moves  it  slightly  southwestward).  The  ring’s  shape  remains 
persistently  circular  and  the  eastward  position  of  the  Loop 
Current  is  misplaced.  NLOM  1/16°  does  show  a  reattachment 
of  the  ring,  but  the  position  of  both  the  ring  and  Loop  Cur¬ 
rent  do  not  agree  with  the  chlorophyll  measurements.  Again, 
NLOM  1/32°  shows  the  best  agreement  with  the  SeaWiFS 
imagery,  with  the  eastward  edge  of  the  eddy  clearly  delin¬ 
eated  by  a  long  plume  of  high  chlorophyll  originating  from 
the  Louisiana  coast.  NLOM  1/32°,  however,  fails  to  place 
the  eastward  frontal  position  of  the  Loop  Current  properly 
between  Cuba  and  the  Florida  Keys.  HYCOM  1/12°  is  the 
only  system  that  accurately  positions  that  specific  front. 
HYCOM  1/12“  also  places  the  reattached  eddy  reasonably 
well,  but  is  missing  some  of  its  lateral  extent,  especially 
the  eastern  area  next  to  the  high  chlorophyll  plume.  NLOM 
1/32°  and  HYCOM  1/12“  continue  to  represent  the  cyclonic 
eddy  depicted  by  high  chlorophyll  on  the  south  side  of  the 
Loop  Current  ring.  This  area  of  high  chlorophyll  is  overlain 
by  the  NCOM  1/8°  Loop  Current  ring  and  is  bisected  by  the 
southern  side  of  the  NLOM  1/16°  Loop  Current  ring. 

08  August  2003:  With  most  of  the  ring  reattached,  the  Loop 
Current  is  now  elongated  and  extends  farther  to  the  west  and 
to  the  north  (Plate  4).  A  fragment  of  the  Loop  Current  ring 
remains  detached  west  of  the  reattached  portion  (small  dark 
area),  but  is  not  depicted  by  any  of  the  systems.  In  NLOM 
1/16",  the  ring  did  not  remain  attached  and  moved  westward 
of  90“W.  The  Loop  Current  position  is  reasonably  well  repre¬ 
sented  east  of  88“W.  The  NCOM  1/8°  SSH  field  shows  very 
little  agreement  between  the  chlorophyll  measurements  and 
the  circular  ring  formed  on  09  July  2003,  now  located  farther 
to  the  west.  In  addition,  the  Loop  Current  does  not  extend  far 
enough  to  the  north.  Both  HYCOM  1/12"  and  NLOM  1/32” 
do  a  good  job  at  capturing  the  full  northwestward  extent  of 
the  Loop  Current.  There  are  some  small  differences  in  the 
representation  of  the  recaptured  Loop  Current  ring.  The  ring 
in  HYCOM  1/12“  is  broader  and  shows  more  closed  SSH  con¬ 
tours,  but  both  have  nearly  the  same  the  ring  center  location. 
Both  HYCOM  1/12“  and  NLOM  1/32“  fail  to  correctly  place 
the  eastward  frontal  position  of  the  Loop  Current,  which  is 
well  delineated  by  high  chlorophyll  in  the  observations. 
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Plate  3.  SSH  from  four  ocean  data  assimilative  systems  on  26  July  2003  (contour  interval  of  10  cm  and  normalized 
over  the  displayed  area),  overlying  ocean  color  imagery  from  SeaWiFS  (natural  log  of  the  concentration  in  mg  nr1). 
The  model  land-sea  boundary  is  defined  by  the  black  contour  line. 
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5.  DISCUSSION 

The  comparisons  performed  in  section  4  illustrate  the 
value  of  the  SeaWiFS  imagery  in  assessing  the  ability  of 
three  ocean  data  assimilative  systems  to  map  mesoscale 
ocean  features  in  the  Gulf  of  Mexico.  For  routine  evaluation, 
quantitative  measures  of  the  frontal  position  of  the  Loop 
Current  could  be  established,  similar  to  what  is  presently 
done  routinely  with  satellite  SST  measurements.  Specifically, 
these  SeaWiFS  measurements  can  provide  additional  infor¬ 
mation  on  the  frontal  position  of  the  Loop  Current,  especially 
during  the  summer  months  when  the  surface  gradient  in  SST 
is  weak  or  nonexistent. 

It  is  apparent  from  the  comparison  that  horizontal  resolu¬ 
tion  is  of  importance  in  properly  representing  such  main 
features  of  the  Gulf  of  Mexico  as  the  Loop  Current  and 
associated  rings,  since  the  NCOM  1/8°  system  does  not 
perform  as  well  as  the  NLOM  1/16°  system  despite  assimi¬ 
lating  the  NLOM  1/16°  SSH  fields.  The  NCOM  1/8°  system 
is  clearly  outperformed  by  the  other  three  systems,  with  the 
high  resolution  NLOM  1/32"  exhibiting  the  highest  skill  and 
the  HYCOM  1/12°  system  running  a  close  second.  HYCOM 
1/12°  performs  better  overall  than  NLOM  1/16°  despite  the 
fact  that  they  have  nearly  the  same  horizontal  resolution  (~8 
km)  in  the  Gulf  of  Mexico. 

In  order  to  illustrate  the  impact  of  horizontal  resolution  and 
of  assimilating  different  SSH  fields  in  NCOM,  we  now  dis¬ 
cuss  the  results  (Plate  5)  of  the  NCOM  system  using  a  1/24" 
Intra-Americas  Sea  configuration  (~  6  km  grid  resolution) 
[Ko  et  al.,  2003], embedded  in  the  NCOM  1/8°  system.  The 
NCOM  1/24°  system  uses  the  same  data  assimilation  tech¬ 
nology  as  NCOM  1/8°,  except  that  it  assimilates  the  model 
independent  MODAS  2-D  SSH  analyses  [Fox  et  al.,  2002] 
instead  of  the  SSH  fields  from  the  NLOM  1/16°  system. 
Like  the  other  systems,  the  forcing  is  from  NOGAPS.  The 
open  boundary  conditions  are  provided  by  NCOM  1/8°  (one¬ 
way  coupling)  and  include  sea  surface  elevation,  transport, 
temperature,  salinity,  and  velocities.  Compared  to  NCOM 
1/8°,  NCOM  1/24°  shows  an  improved  representation  of  the 
Loop  Current  ring  shedding  event  on  all  four  days  (Plate  5): 
the  separated  ring  does  not  remain  as  circular  as  it  was  the 
case  in  NCOM  1/8°  (although  it  is  still  more  circular  than 
depicted  by  SeaWiFS)  and  it  is  in  generally  at  the  right  loca¬ 
tion  at  the  end  of  the  sequence  (though  too  far  south  by  half  a 
degree).  However,  at  this  time,  NCOM  1/24°  has  insufficient 
northward  penetration  of  the  Loop  Current  (less  even  than 
NCOM  1/8°). 

Overall,  Plates  1-5  demonstrate  striking  differences  in 
the  mapping  of  the  Loop  Current  and  associated  eddies.  One 
could  have  expected  a  closer  agreement  between  the  systems 
given  the  fact  that  they  all  assimilate  SSH  fields  from  the 


same  altimeters.  In  order  to  further  assess  the  impact  of  the 
data  assimilation  techniques  used  by  the  systems  as  well 
as  the  ability  of  ocean  color  to  track  the  observed  SSH  fea¬ 
tures,  model-independent  SSH  analyses  (except  for  the  SSH 
mean)  performed  by  the  Colorado  Center  for  Astrodynamics 
Research  at  the  University  of  Colorado,  Boulder  [ Leben 
et  al.,  2002]  are  presented  in  Plate  6  for  comparison.  The 
maps  are  produced  from  Topex-Poseidon,  JASON,  GFO,  and 
ENVISAT  altimeter  data  anomalies  added  to  a  mean  SSH 
field  based  on  a  CUPOM  model  simulation  from  1993-1999 
(see  chapter  by  Kantha  et  al.  in  this  volume  for  details  on 
CUPOM).  These  hindcast  analyses  have  the  potential  for 
higher  accuracy  because  they  include  altimeter  data  from 
two  additional  satellites,  Topex-Poseidon  and  ENVISAT, 
not  then  available  in  the  real-time  data  stream  at  the  Naval 
Oceanographic  Office,  the  data  stream  used  by  all  five  of 
the  data-assimilative  models  compared  with  the  ocean  color 
imagery.  Although  less  precise  without  the  overlay,  com¬ 
parison  of  Plate  6  with  the  ocean  color  imagery  of  Plates  1-5 
clearly  shows  that  ocean  color  and  the  SSH  analyses  track 
the  same  major  ocean  features  present  in  the  Gulf  of  Mexico 
in  quite  good  agreement. 

Horizontal  resolution  of  the  ocean  model  is  a  factor  in  the 
accuracy  of  the  mapping  by  the  data  assimilative  systems,  but 
HYCOM  1/12°  also  outperforms  the  NLOM  1/16°  and  NCOM 
1/24°,  which  have  similar  or  finer  resolution,  showing  that 
other  factors  (some  influenced  by  resolution)  are  important 
as  well.  Other  differences  in  the  prediction  systems  include 
the  vertical  coordinate,  the  vertical  resolution,  the  bottom 
topography,  and  the  data  assimilation.  For  data  assimilation, 
the  assimilation  of  SSH  and  the  downward  projection  of  the 
SSH  updates  are  particularly  significant.  The  NLOM  systems 
assimilate  altimeter  track  data  and  use  the  statistical  inference 
technique  of  Hurlburt  et  al.  [1990]  for  downward  projection 
which,  unlike  the  other  systems,  includes  projection  into  the 
abyssal  ocean,  a  potential  advantage  in  representing  baroclinic 
instability.  HYCOM  1/12°  uses  the  Cooper  and  Haines(1996) 
scheme  for  downward  projection  of  SSH  updates,  while  the 
NCOM  systems  assimilate  SSH  (and  SST)  indirectly  via 
MODAS  synthetic  T  &  S  profiles  [Far  et  al.,  2002].  NCOM 
1/8”  assimilates  the  NLOM  1/16°  SSH  fields  while  HYCOM 
1/12°  and  NCOM  1/24°  assimilate  the  model-independent 
MODAS2D  analyses  of  SSH  anomalies  [Fax  et  al..  2002] 
operational  at  NAVOCEANO,  but  use  a  different  mean  SSH,  a 
MICOM  mean  from  an  atmospherically-forced  simulation  for 
HYCOM  and  the  surface  dynamic  height  from  the  MODAS 
hydrographic  climatology  for  NCOM  1/24°. 

The  simulation  skill  of  the  ocean  model  without  data  assimi¬ 
lation  also  strongly  influences  prediction  system  skill.  This 
skill  is  influenced  by  model  resolution  and  other  factors, 
including  those  listed  earlier.  In  the  Gulf  of  Mexico,  accurate 
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Plate  5.  SSH  from  the  1/24"  NCOM  data  assimilative  system  oil  09.  19.  26  July,  and  08  August  2003  (contour  interval 
of  10  cm  and  normalized  over  the  displayed  area),  overlying  ocean  color  imagery  from  SeaVViFS  (natural  log  of  the 
concentration  in  mg  nr  ').  The  model  land-sea  boundary  is  defined  by  the  black  contour  line. 
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inflow  through  the  Yucatan  Channel  is  critical  and  may  not  be 
strongly  constrained  by  the  data  assimilation.  Accurate  inflow 
transport  and  realistic  vertical  and  horizontal  structure  of  the 
inflow  are  needed,  as  this  influences  the  northward  Loop 
Current  penetration  into  the  basin  and  the  stability  of  the  Loop 
Current  after  separation  from  the  Campeche  Bank.  The  inflow 
includes  strong  surface  trapping  with  a  high  speed  core  on  the 
west  side  of  the  channel.  The  current  structure  in  HYCOM 
1/12°  [Chassignet  et  al„  2005]  has  shown  excellent  compari¬ 
son  with  observations  [Sheinbaum  et  al.,  2002;  Abascal  et  al., 
2003]  and  a  strong  surface  jet  on  the  west  side  of  the  channel 
is  indicated  by  the  tightly  packed  SSH  contours  in  the  figures. 
NLOM  1/16°  and  NLOM  1/32°  show  a  similar  SSH  signature 
for  the  inflow.  NLOM  simulations  without  data  assimilation 
give  higher  inflow  transport  and  greater  surface  trapping 
at  1/32°  resolution.  The  NCOM  systems  have  lower  inflow 
transport,  an  indication  of  why  the  HYCOM  and  NLOM 
systems  tend  to  give  better  agreement  with  SeaWiFS  for  the 
northward  Loop  Current  penetration.  Excessive  abyssal  inflow 
on  the  west  side  of  the  channel  can  also  inhibit  Loop  Current 
penetration  [. Hurlburtand  Thompson,  1980].  NLOM  1/16°  has 
the  opposite  problem,  with  an  unrealistic  westward  abyssal 
current  along  the  northern  continental  slope  due  to  south¬ 
ward  leakage  of  some  flow  from  the  Atlantic  Deep  Western 
Boundary  Current  (DWBC)  through  the  Florida  Strait  (which 
does  not  occur  in  NLOM  1/32°  after  modification  of  the  bot¬ 
tom  boundary  condition).  Abyssal  currents  can  advect  upper 
ocean  currents  where  they  intersect  at  large  angles  [Hurlburt 
and  Thompson,  1980;  Hurlburt  et  al,  1996].  This  causes  the 
Loop  Current  ring  to  move  westward  too  rapidly  (Plate  4). 

Another  key  difference  in  the  systems  is  their  ability  to 
represent  baroclinic  instability,  a  phenomenon  manifested  by 
contortions  of  the  Loop  Current  and  associated  eddies  (e.g., 
noncircular  eddies)  and  by  cyclonic  eddies  on  the  periphery 
of  the  Loop  Current  and  associated  anticyclonic  rings.  This 
is  particularly  well  represented  in  the  NLOM  1/32°  system 
and  is  also  exhibited  in  the  NLOM  1/16°  system.  It  is  seen 
in  the  HYCOM  1/12°  system,  but  appears  to  be  inhibited  by 
the  assimilation  of  the  relatively  smooth  model-independent 
SSH  anomaly  analyses.  The  HYCOM  1/12°  and  NLOM  1/32° 
systems  also  depict  the  cyclonic  eddy  on  the  south  side  of  the 
Loop  Current  ring  (Plates  2-3)  which  is  missing  in  the  NCOM 
1/8°  and  NCOM  1/24°  systems  and  is  poorly  represented  in 
the  NLOM  1/16°  system.  In  the  two  NCOM  systems,  the  Loop 
Current  eddy  is  too  circular  and  the  observed  contortions  of 
the  Loop  Current  and  associated  eddies  are  limited,  especially 
in  NCOM  1/8°  where  the  Loop  Current  ring  is  much  more 
circular  than  the  NLOM  1/16°  ring  it  assimilated  (Plates  1-4). 
Particularly  striking  are  the  differences  in  HYCOM  1/12° 
and  NCOM  1/24°,  which  both  assimilated  the  same  model- 
independent  SSH  anomaly  analyses.  These  are  associated 


with  the  differences  in  model  dynamics,  choice  of  mean  SSH, 
and  the  related  Yucatan  inflow  (see  chapter  by  Oey  et  al.  in 
this  volume  for  a  comparison  between  the  observed  Yucatan 
Channel  flow  structure  and  its  representation  in  six  different 
models).  The  HYCOM  1/12°  results  are  actually  more  simi¬ 
lar  to  the  NLOM  1/32°,  the  two  in  best  agreement  with  the 
SeaWiFS  imagery. 

6.  CONCLUSION 

The  results  of  this  study  illustrate  the  value  of  SeaWiFS 
ocean  color  imagery  in  assessing  the  ability  of  five  ocean 
data  assimilative  systems  to  map  mesoscale  variability  in 
the  Gulf  of  Mexico.  A  quantitative  assessment  was  not  per¬ 
formed  as  it  would  have  required  evaluating  the  frontal 
position  of  the  Loop  Current  and  eddies  from  the  SeaWiFS 
imagery,  a  task  beyond  the  scope  of  this  paper. 

The  SeaWiFS  imagery  is  very  effective  in  monitoring  the 
changes  in  the  mesoscale  activity  in  the  Gulf  of  Mexico,  in  dif¬ 
ferentiating  the  ability  of  the  five  systems  to  map  the  evolution 
of  the  Loop  Current  and  associated  warm  and  cold  eddies,  and 
in  helping  to  diagnose  specific  strengths  and  weaknesses  of 
the  systems.  This  was  made  possible  by  improved  compositing 
methods  which  effectively  removed  the  clouds  and  retained 
ocean  features  in  the  surface  chlorophyll  properties.  The 
imagery  provides  a  unique  independent  data  set  for  validating 
numerical  models  by  determining  the  locations  of  mesoscale 
features.  Both  chlorophyll-rich  and  chlorophyll-poor  regions 
proved  useful.  High  chlorophyll  plumes  from  the  Mississppi 
River  were  advected  offshore  and  helped  outline  the  eastern 
edge  of  the  Loop  Current.  Elevated  chlorophyll  was  also  pres¬ 
ent  in  the  center  of  some  cyclonic  eddies.  Dark  chlorophyll- 
poor  areas  clearly  represented  the  interior  of  the  Loop  Current 
and  Loop  Current  rings. 

In  addition,  the  study  clearly  illustrates  that  biological 
responses  of  the  surface  waters  are  strongly  linked  to  the 
physical  events  and  processes.  It  also  shows  the  increasing 
ability  of  eddy-resolving,  data  assimilative  ocean  models  to 
map  ocean  features  and  events  that  are  biologically  relevant. 
This  capability  can  help  to  link  observed  biological  evolution 
quantitatively  to  observed  circulation  processes  and  events 
such  as  upwelling,  subduction,  advection,  convergence, 
divergence  and  vertical  mixing,  a  topic  that  should  be  the 
subject  of  future  research. 
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